Level density of a Fermi gas: average growth and fluctuations 
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We compute the level density of a two-component Fermi gas as a function of the number of 
particles, angular momentum and excitation energy. The result includes smooth low-energy cor- 
rections to the leading Bethe term (connected to a generalization of the partition problem and 
Hardy-Ramanujan formula) plus oscillatory corrections that describe shell effects. When applied 
to nuclear level densities, the theory provides a unified formulation valid from low-lying states up 
to levels entering the continuum. The comparison with experimental data from neutron resonances 
gives excellent results. 

PACS numbers: 03.75.Ss, 21.10.Ma, 24.60.-k 



Many physical properties of an interacting Fermi gas 
depend on the number of available states at a given 
energy, like for instance the optical and electrical re- 
sponse of solids, or the reaction rates in nuclear pro- 
cesses. The most common framework to compute the 
many-body (MB) density of states (DOS) is a mean- 
field approximation, where each (quasi-)particle moves 
independently in an average self-consistent potential. In 
this case, the energy of the Fermi gas is expressed as the 
sum of the occupied single-particle (SP) energies. The 
computation of the MB DOS is thus reduced to a combi- 
natorial problem: to count the different ways into which 
the energy can be distributed among the particles. A 
first answer to this question was given by Bethe [2, who 
showed that at high excitation energies Q (compared to 
the SP spacing at Fermi energy ep) and for two types of 
fermions (protons and neutrons), the MB density grows 
like /9mb(Q) exp {2y/aQ^ /Q^^^. The Fermi gas param- 
eter a = 7r^p(eF)/6 depends only on the average SP DOS 
at ep, p(eF). 

In practice the parameter a is often used as a fitting 
parameter. For a given excitation energy Q and parti- 
cle number A, a{Q,A) is extracted from the available 
experimental data. In this way, important deviations 
from the independent particle model predictions are ob- 
served. Though there are certainly effects that are be- 
yond that model, our purpose here is to show that a 
detailed treatment is able to describe features of indi- 
vidual systems with good accuracy, therefore providing 
a solid theoretical basis for extrapolations to unknown 
regions and for improvements. Generalizing the results 
obtained in Ref . ,2| to a two-component system of given 
angular momentum, we show that Bethe's result can be 
viewed as the first (smooth) term of an expansion. The 
corrections to that term do not enter as corrections of the 
a-parameter (as was often assumed in the past), but sim- 
ply as additional terms in the exponential. A first series 
of terms are smooth in A and Q, and provide higher- 
order (in inverse powers of the excitation energy) cor- 
rections. Keeping only the first correction generates a 
uniform expression, which cancel the divergence produce 



by the Q"^/* term at low energies, and therefore make 
unnecessary the use of composite models (a la Gilbert- 
Cameron). On top of the smooth contributions are os- 
cillatory terms, that describe density fluctuations as A 
varies. These are shell effects, which turn out to be re- 
lated to the fluctuations of the total energy of the system. 
A detailed description of these fluctuations and of the 
relevant energy scales is provided. Finally, a comparison 
of the results with the nuclear level density at neutron 
threshold is made. With a few adjustable parameters, a 
very good overall agreement is obtained, with a relative 
error ~ 10% for the logarithm of the density of the 295 
nuclei analyzed. 

From a theoretical point of view, Ref. ^2] and the 
present work may be viewed as a generalization to MB 
systems of the theory developed to describe the SP DOS 
0, S ■ Compared to direct generalizations of the SP 
Weyl's law and of semiclassical trace formulas 0, the 
statistical approach employed here has the advantage of 
leading to simple and physically meaningful results. 

The DOS at energy E oi a. system composed of Z pro- 
tons, N neutrons and with projection AI of the angular 
momentum on some given axis is defined as 

p„^{E, N, Z, M) ^Y,5{E~ E'')5{N - AJ^) 

5{Z - A^'z) 5{M ~ M") . (1) 

The index v denotes all the possible neutron and pro- 
ton SP configurations (of arbitrary number of particles) , 
yl^ — ^ are the neutron and proton number of par- 

ticles, respectively (n^ ^ = 0, 1 are the corresponding oc- 
cupation numbers of the i-th SP state, and A = or 
Py.E'^^ Ea and = are 

the energy and angular momentum projection, where e\^i 
and m\^i denote the SP energies and angular momentum 
projections, respectively. 

The conservation of the angular momentum projection 
is introduced in order to deal with only a subset of states, 
those with given total angular momentum J. A standard 
treatment 0,111 of this degree of freedom leads to a DOS 
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of angular momentum J of the form 



p«u{Q,N, Z, J) 



2J+1 + 



2v/2^cr3 



Pmb((3, n, z), 



(2) 

where p^i^{Q,N,Z) is the total MB DOS, Q = E - Eq 
is the excitation energy measured with respect to the 
ground state of the system, and a is the spin cutoff pa- 
rameter. 

For an arbitrary SP spectrum the computation of the 
density of excited states is a difficult combinatorial prob- 
lem for which no exact solution exists. There is, how- 
ever, a particular case that can be worked out explicitly: 
when the SP spectrum consists of equidistant levels sep- 
arated by 5 (we assume, for simplicity, that the neutron 
and proton spacings are equal). The MB excitation en- 
ergies are then given by the sum of two integers corre- 
sponding to the total energy of each of the components, 
Q — {j + k) S — K S. Each MB state characterized by 
an integer K has a nontrivial degeneracy. The compu- 
tation of the degeneracy reduces to the computation of 
the number of ways into which the total energy may be 
distributed among the two components, and of the dif- 
ferent ways the partial energy of each component can be 
distributed among its elements. This leads to compute 
the value of the function P2{K) — J2f=iPij) ~ j)' 
where p{j) is the partition of j (the number of ways into 
which the integer j can be decomposed as a sum of inte- 
gers). We are assuming here, to avoid finite size effects, 
that the excitation energy is small compared to the Fermi 
energy of each component. Based on the work of Hardy 
and Ramanujan, an exact expression (written as a con- 
vergent series) for p{j) was obtained by Rademacher 0. 
We have adapted their method (i.e., the circle method, cf 
Ref . [lOj ) to obtain an exact formula for p2{K). Putting 
back the appropriate units, the MB density can be ex- 
pressed in terms of p2 as Pmb/p = 2^/^p2{K = p (9/2), 
where p = pp -I- pjy = 2/5 is the total (proton -I- neu- 
trons) SP average density. Then, expressing the exact 
result as an expansion in terms of p Q valid in the range 
< Q < N5, Z5, we obtain 



,(Q,iV,Z)/p = 



61/4 



12 (p g)'V4 



(3) 



where the "entropy" S = Seq of the equidistant spectrum 
is given by 




Se. 



plus 0((p Q)^^^"^) corrections that can be computed but 
are not given here. The prefactor in Eq.Q and the first 
term of the expansion |0J reproduce Bethe's formula Ij. 
The additional terms provide further smooth corrections 
of higher order in inverse powers of the excitation energy. 
Though Eq.Q represents an asymptotic expansion, we 
find that an accurate uniform approximation is obtained 
by keeping only up to the l/^/pQ term that kills, when 



Q ^ 0, the divergence produced by the (pQ)~^^^ in the 
prefactor. It is also interesting to note that the correct 
coefficients of the correction terms in Q are obtained 
through the expansion of the exact result, whereas a sad- 
dle point approximation of the sum involved in p2{K) 
leads to wrong coefficients. 

The previous expression describes in detail the MB 
DOS for a SP spectrum made of equidistant levels. How- 
ever, it is clearly unphysical in most situations. A generic 
SP spectrum contains fluctuations, which are manifested 
at the scale of the average distance between levels, but 
also on much larger scales (see Ref.[ll| for a review). 
What is missing in Eq. Q are the fluctuations in the MB 
density induced by the SP fluctuations. In this respect, 
one may consider Eq.Q as the MB analog of the Weyl 
or Wigner-Kirkwood expansions. 

It remains to compute the MB level density for an ar- 
bitrary SP spectrum, including fluctuations. The way to 
do it was shown, for a single-component gas, in Ref.0]. 
The method uses a saddle point approximation of the 
inverse Laplace transform of the MB density. We have 
adapted that calculation, following similar lines, to treat 
Eq.Q, that includes two components and angular mo- 
mentum conservation. The result may be written under 
the form of Eqs.|21) ^^i^ (jSJ, but with the entropy in the 
latter equation given by 



S — Sea + 7^ 



£{N, Z,0)-£{N, Z, T) 



(5) 



The parameter T is the temperature, connected to the 
excitation energy Q through the usual relation Q = 

^2pTV6 = a T\ £{N, Z, T) = Ea / de M^) ^ fi^, T) 
is the fluctuating part of the energy of the system at 
temperature T and chemical potential p, ~ ep fixed, ne- 
glecting temperature variations, by the particle-number 
conditions N ^ Z = de p\{e). The function px{e) = 
J2j '5(e — is the SP density of the component A, and 

PA(e) = Px{^) ~Pa(^) its fluctuating part. £{N,Z,0) is 
thus the fluctuating part of the ground state energy of 
the system. Finally, Seq in Eq.© is given by Eq.Q), 
with p the total average SP density of the system at 
Fermi energy. In fact, for an arbitrary spectrum the sad- 
dle point technique does not allow to derive the terms 
of order (pQ)^^^^ and higher in Seq- The corrections ob- 
tained from an equidistant spectrum are thus conjectured 
to provide a good approximation to the corrections of 
the smooth part of an arbitrary system, but the validity 
of this statement has to be conflrmed. An explicit nu- 
merical verification of its validity for a two-dimensional 
one-component system was done in Ref. . 

The function £{N,Z,T) presents oscillations when N 
or Z are varied, in contrast to the more gentle variations 
as a function of T (a detailed description of the fluctua- 
tions and of their T-dependence is given below) . The MB 
level density contains now two types of terms: some that 
vary smoothly, and others that fluctuate as the number of 
particles varies. The result presented above contains the 
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dominant smooth and oscillatory terms. In the deriva- 
tion of Eqs.Q, ® and © we have neglected other terms 
(for instance, the chemical potentials and T have small 
corrections that depend on J, and thus strictly speak- 
ing the factorization |(5J) of the angular momentum is not 
exact, etc). A detailed account of the derivation will be 
given elsewhere. 

It is remarkable that the MB level density at excita- 
tion energy Q depends explicitly on the ground-state en- 
ergy fluctuations £{N,Z,0). A convenient way to ana- 
lyze the behavior of the fluctuating part of the entropy 

S = £{N, Z, 0) - £{N, Z, T) /T is through a semiclas- 

sical theory. The result is an expression for S written as 
a sum over all the classical periodic orbits of the mean 
field potential. The main conclusions that can be drawn 
from that expression are now listed. To be specific, we 
consider the particular case of an atomic nucleus of Z pro- 
tons and N neutrons: (a) as the mass number A — Z + N 
varies at a fixed excitation energy Q, S presents oscilla- 
tions of characteristic period 5A « (7r/3)A^/'^, that are 
independent of Q; (b) when Q varies at fixed particle 
number A, S does not present similar oscillations, but 
rather gentle variations; (c) the typical amplitude <Tg of 
S at given [Q, A) depends on the dynamical properties of 
the classical dynamics (integrable or chaotic); (d) since 
the mean field dynamics of most nuclei is well approx- 
imated by a regular motion then the behavior of 
the typical amplitude of S is given, to a first approxima- 
tion, by those of a regular dynamics, that we now detail; 
(e) using the definition of the temperature T = yjQ/a 
with a = Tp--pl% K, AjXh MeV^^ |i3|, we find that there 
is only one relevant temperature scale in the variation 
of (j-g with T, given by w 4/^1/3 MeV; for conve- 
nience we also introduce = (27r^p)^^ ~ 1.3/A MeV 
(the temperature associated to the SP mean level spac- 
ing), and g — Tc/Ts ~ 3^^/"^; (f) at low temperatures, 
~ ^JT/Ts\ the typical amplitude of shell effects in 
the MB DOS therefore increases from at T = up to 
~ y/g at T ~ Tc] (g) at temperatures of order the 
amplitude is maximal, and starts to decrease for T > Tc] 
(h) in the limit T ^ the typical amplitude tends to 
zero as « a/ g/6 Tc/T; using the previous values of 
Tc and g this gives ~ 2\/2/T MeV; we thus predict a 
slow power-law decay of the amplitude of shell effects at 
high temperatures. 

We now turn to a direct application of the previous re- 
sults to experimental data. Though it will be important 
to make a systematic analysis of the validity of Eqs.Q, 
Q and Eq.lO) and of their predictions at different en- 
ergies and mass numbers, we restrict here to a compar- 
ison with slow neutron resonances, which have been ex- 
perimentally studied for a large number of nuclei [l3l| . 
The excitation energies of neutron resonances coincide 
with the neutron binding energies, Q = Qn = Sn{N, Z), 
whose values are in the range 6 — 8 MeV for most nu- 
clei. This corresponds to a temperature T„ « 8/^/A 



MeV. According to the previous results, the typical am- 
plitude of the fluctuations depends on temperature, with 
a maximum at T « Tc. At neutron resonances the ratio 
Tn/Tc « 2/^1/6. From A = 30 to A ^ 250, this ratio 
varies from 1.13 to 0.8. We thus find that at excitation 
energies Q ~ Qm the temperature is very close to T^, 
shell effects are maximal. We expect a typical value of 
S{Q,N,Z) very close to its maximum ~ y/g « ^/SA^^"^ 
(this varies from 5.4 to 11 in the previous range of A). In 
contrast, in the same particle-number range the first cor- 
recting term (proportional to (p Q)~^^^) in the smooth 
expansion 10} varies from 0.32 to 0.11. That term, and 
the following ones in the expansion, can thus be neglected 
at Q ~ Qn- 

To make a comparison with experiments we need the 
different quantities involved in the theoretical expres- 
sions. One possibility is to compute them from a particu- 
lar model. In our case, however, in order to avoid model- 
dependent features and to make a direct test of our pre- 
dictions we prefer to extract as much information as pos- 
sible from experimental data. For each nucleus, the exci- 
tation energy at neutron threshold Q = Qn = Sn{N, Z) 
is taken from the experimental value of Sn{N,Z), and 
Tn{N, Z) = ^jQn/a. For £{N, Z, T„) we have no exper- 
imental data available. We compute it as follows. Semi- 
classically, £{N, Z, Tn) is written as a sum over the peri- 
odic orbits p (and repetitions) of the mean field potential 
[lll | . The analysis of the temperature dependence of that 
sum and of the main contributing orbits leads to the ap- 
proximation £{N, Z, Tn) ~ K„ £{N, Z, 0), where k„ is the 
average over the shortest periodic orbits p of the function 
K{xp) = Xp/ sinh(xp), where Xp = iiripA^/^Tn/ {'ie-p) and 
ip is the length of the periodic orbit p measured in units 
of the nuclear radius (notice the mass number and tem- 
perature dependence of Xp). For each nucleus, A and 
Tn are given and the average k„ = K(iV, Z, T„) is com- 
puted. In practice, the average is estimated using the 
shortest periodic orbits of a spherical cavity of radius 
R = 1.2Ai/3fm. 

The expression of the entropy S in Eq.|j2Jl takes now 
the form 

S{Qn,N, Z) = 2v/^ + (1 - ^) £{N, Z, 0)/r„ . (6) 

Finally, £{N,Z^Q) and a are required. £(iV, Z, 0) can 
be obtained from the experimental value of the ground- 
state energy. It corresponds to the fluctuating part of the 
nuclear binding energy, that we compute by subtracting 
from the 1995 Audi-Wapstra compilation Ti'l the liq- 
uid drop expression £ — a.t,A — agA^^^ — QcZ"^ /A^/'^ — 
aA{N — Z)'^/A, using the parameters (from Ref. 15] ) 
Qy = 15.67, as = 17.23, ua = 23.29, and = 0.714 (all 
in MeV; we have moreover excluded the pairing term). 
This parametrization produces a fluctuating part whose 
average (over A) {£{N, Z,0))a ~ 0. However, the de- 
termination of the average of the fluctuating part is a 
delicate question that deserves a careful discussion. Due 
to the discrete variation of the chemical potential as 
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the mass number varies, one can verify that generically 
(£(-/V, Z, 0))^ is non-zero. We have therefore added to the 
fluctuating part a term bA + c, where b and c are two con- 
stants. Equation (j^ thus depends on three constants, a, 
&, and c, that we fix by minimizing the root mean square 
error with respect to the experimental value of the den- 
sity, Sexp (obtained by computing S from Eqs.® and 
Q when Pmb{Q, N, Z, J) is the experimental DOS, J the 
ground state angular momentum and cr^ « O.lbaA^^'^Tn 
The result is a = A/WA2 MeY'^, b = -0.019 
MeV and c = 7.9 MeV. The comparison is made in Fig.l. 
The experimental values Sexp shown on the top part are 
to be compared with the "theoretical" entropies plotted 
in the middle part. A clear overall agreement is observed. 
For most nuclei, the relative error in the lower panel is 
smaller than 10%, with some remaining structure as a 
function of A, and larger deviations for closed shells (we 
suspect that this is due to our very schematic estimate 

of K„). 

The precision of the present calculation, with only 
three adjusted parameters, is comparable to the best fits 
obtained nowadays. We can in fact make the compari- 
son more precise by noticing that Eq.© can be approx- 
imated, using an effective value of a, by S{Qn, N, Z) « 

2y/aef Qn, where a^f ^ a [1 + £{N,Z, 0){1 - k)/Q„]. 
Under this form, Eq.® is quite similar to one of the 
best phenomenological formulas studied so far, proposed 
by Ignatyuk and collaborators [TsL . 

To conclude, we have derived an explicit formula for 
the MB DOS of a two component Fermi gas of fixed an- 
gular momentum. The results were applied to the partic- 
ular case of nuclear level densities, where precise predic- 
tions for the smooth dependence and shell fluctuations 
as a function of excitation energy and mass number were 
made. Good agreement between theory and experiment 
in the region of neutron resonances is found. Although 
it was derived within an independent particle model, the 



comparison with experiments shows that the flnal result 
is probably of more general validity and includes, through 
the energy fluctuations, effects like pairing. Going to high 
excitation energies, the main prediction is the decay of 
shell effects when Q ~ Qn (with a power-law tail) . How- 
ever, that prediction is valid for closed systems. Before 
proceeding in that direction, the theory should be im- 
proved to include finite size effects (e.g. a finite number 
of nucleons) as well as the influence of the continuum. 
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FIG. 1: Entropy <S as a function of the mass number A for 
nuclear level densities at neutron threshold, (a) experimental 
values; (b) theoretical prediction; (c) relative error. 
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